G^ 



00 



Oh! 



Integrable Hierarchies and Contact Terms in u-plane 
Integrals of Topologically Twisted Supersymmetric 

Gauge Theories 



Kanehisa Takasaki 



X5 

D ■ Department of Fundamental Sciences, Kyoto University 

OO ■ Yoshida, Sakyo-ku, Kyoto 606, Japan 

El-mail: tcLkasaki@yukawa.kyoto-u.ac.jp 

^ ' Abstract 

. The w-plane integrals of topologically twisted N = 2 supersymmetric gauge theo- 



ries generally contain contact terms of nonlocal topological observables. This paper 
proposes an interpretation of these contact terms from the point of view of inte- 
grable hierarchies and their Whitham deformations. This is inspired by Mariiio and 
Moore's remark that the blowup formula of the n-plane integral contains a piece 
that can be interpreted as a single-time tan function of an integrable hierarchy. This 
^ ■ single-time tan function can be extended to a multi-time version without spoiling 

Cd I the modular invariance of the blowup formula. The multi-time tau function is com- 

prised of a Gaussian factor e*^*^*^'*^' ' ) and a theta function. The time variables t„ 
play the role of physical coupling constants of 2-observables In{B) carried by the 
exceptional divisor B. The coefficients qmn of the Gaussian part are identified to be 
the contact terms of these 2-observables. This identification is further examined in 
the language of Whitham equations. All relevant quantities are written in the form 
of derivatives of the prepotential. 
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1 Introduction 



The "u-plane integral" of Moore and Witten gives an exact answer to the contribution 
of the Coulomb branch in correlation functions of topologically twisted four- dimensional 
N = 2 supersymmetric gauge theories. Moore and Witten considered the case of SU{2) 
and 5*0(3) only. Marino and Moore extended the w-plain integral to more general, 
higher rank gauge groups. Losev et al. studied the w-plane integral from a somewhat 
different point of view ( "Gromov- Witten paradigm"). 

The correlation functions of those topologically twisted N = 2 supersymmetric gauge 
theories are the Donaldson- Witten invariants of the four- manifold X . Their generating 
function is given by a path integral of the form 



with book-keeping variables x and y coupled to the 2-cycle 5* G H2{X,Z) and the 0- 
cycle O G Hq[X,Z). I{S) and 0{P) are associated observables ( "2-observable" and 
"0-observable" ) . The Coulomb branch has a nonvanishing contribution only if < 1, 
and this contribution is given by an integral of the following form ("u-plane integral") 
over the Coulomb moduli space: 



X and a are respectively the Euler characteristic and the signature of X; f/ is a contribu- 
tion of the 0-observable; T is the " contact term" induced by the 2-observable; \E' is the 
photon partition function, which is a lattice sum of the Siegel-Narain type over the tensor 
product of H2{X, Z) and the weight lattice of the gauge group G. Since the geometry of 
the Coulomb moduli space is determined by the low energy effective action of Seiberg and 
Witten ( "Seiberg- Witten solution") [^, this w-plane integral gives a complete answer to 
the physics of the Coulomb branch at least in the topologically twisted theories. 

The low energy effective action of the Coulomb branch is generally related to an 
integrable Hamiltonian system P|. Our subsequent consideration is mostly focussed on 
the SU (N) theory without matter hypermultiplet {Nf = 0) 0. In this case, the integrable 
Hamiltonian system is the A^-periodic Toda chain. Its spectral curve can be written 




(1) 




(2) 



- A-^P{x)z + 1 = 0, 



(3) 
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where P{x) is a polynomial of the form 

N 

P{x)=x'' -Y.Ujx''-'. (4) 

i=2 

This hyperelliptic curve of genus g = N—1 now plays the role of the Seiberg-Witten elliptic 
curves in the SU{2) theories ||^. The coefficients Uj, which give a Poisson-commuting set 
of Hamiltonians of the Toda chain, are thereby identified with the Coulomb moduli. A is 
the energy scale of the low energy theory. The low energy effective action is written in 
terms of a prepotential JF = jF(a), a = (ai, ■ ■ ■ , a]\r_i). This prepotential is determined 
(though as an implicit function) by the functional relation 

d:F 



of the period integrals 
of the meromorphic differential 



aj = (j) rfS'sw, = f '^'S'sw (6) 



dSsw = xdlogz (7) 

along a symplectic basis aj, Pj (j = 1, ■ ■ ■ , — 1) of cycles on the Riemann surface of the 
spectral curve. The differential dSsw obeys the fundamental relation 

dajdSsw\z=const = duj. (8) 

The left hand side means differentiating dSsw while keeping z constant, duj (j = 
basis of holomorphic differentials normalized as 

duj = 6jk. (9) 

This equation connecting dS^w and duj shows a link with the notion of "Whitham equa- 
tions" for adiabatic deformations of algebro-geometric solutions of integrable systems. 

It is accordingly natural to expect a similar relation of the w-plane integrals to some 
integrable systems and associated Whitham equations. Marino and Moore and Losev 
et al. 10, independently, presented several interesting remarks towards that direction. 
The purpose of this paper is to examine their remarks in more detail. 
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2 Blowup formula and contact terms 



One of remarks of Marino and Moore P] is that the "blowup formula" of the u-plane 
integrals contains a factor that can be interpreted as a special tau function of the Toda 
chain (or, more precisely, the Toda lattice hierarchy [^]). 

The blowup formula is concerned with the blowup X of the four-manifold X at a point 
P. Let B denote the homology class of the exceptional divisor (i.e., the inverse image of 
the blowup point), and take the homology class S = S + tB with a parameter t on X. 
Furthermore, let the metric of X be such that the Poincare dual of B is anti-self-dual, 
i.e., -B+ = 0. In the case of G = SU{N) and Nj = 0, the blowup formula shows that the 
effect of blowup is just to replace the 0-observable factor as 



Here a and f3 are some numerical constants, A the discriminant of the above Toda spectral 
curve, and the last factor is the ordinary theta function {g = N — 1) 

Q^^siw \r)= J2 exp [vri < £ + 7, T(£ + 7) > +27ri < £ + 7, w + 5 >] (11) 
eeza 

with the half-characteristic (7, 6) determined by the setup of the w-plane integral. The 
period matrix T = {Tjk)j^k=i,-,g is defined by 

dujj = Tjk- (12) 

13k 

The vector V = (V^)j=i, ■••,£/ is a gradient vector, with respect to a = {ai, ■ ■ ■ , ag), of a 
function V = V{a): 

dV 

This function V appears in the definition of the 2-observable I{S) as 

I{S) = const. / G^P, (14) 



where G is an operator that generates a standard solution of the descent equations for 
observables ^ 0. The contact term T also depends on this potential function V, therefore 
should be written T-p more precisely. 

It is the last two factors in (|T0|) that Marino and Moore identified to be a tau function 
of the Toda lattice hierarchy. This interpretation is very suggestive in the following sense. 
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• The parameter t is interpreted as a time variable in the hierarchy. This is in sharp 
contrast with the status of the Toda chain in the aforementioned description of 
the Coulomb moduli space. The role of the Toda chain is simply to supply a g- 
dimensional family of curves along with a special-geometric structure; the dynamics 
of the Toda chain, as an integrable Hamiltonian system, plays no role. In the it- 
plane integral, meanwhile, the time variable becomes a physical coupling constant 
of the observable I{B). 

• The contact term T and the directional vector V are determined by the potential 
V. In this sense, P is a Hamiltonian of the Toda lattice hierarchy. In fact, any 
polynomial (or holomorphic function) of the Coulomb moduli can be used as V. 
One can consider a set of commuting flows with a set of Hamiltonians Vi,V2, - ■ ■ 
and associated time variables ti,t2," ■ ■• These commuting flows generally form a 
subhierarchy of the standard full Toda lattice hierarchy with two infinite series of 
time variables {n = 1,2, ■ ■ ■) ^ Of course, the notion of tau function is also 
meaningful for such a subhierarchy. This strongly suggests that the blowup formula, 
too, can be extended in that way. 



3 Insertion of more than one 2-observables 

Let us specify the implication of the second point above. The potentials Vi,V2, ■ ■ ■, 
determine the 2-observables 

UB) = const. / G^Vn- (15) 
Jb 

These 2-observables can be used to deform the correlation functions as 

(exp(/(5) + 0(P))) ^ {exp{l{S) + J2tnUB) + 0{P))). (16) 

n 

(The book-keeping parameters are set to a; = y = 1.) This will modify the last two factors 
of (|10D as 

exp{-t^T)Q,4^tV I T) 

^Actually, the Toda lattice hierarchy has another variable to — the lattice coordinates. In order to 
avoid unnecessary complication, we shall not consider it, or just put to — 0. 
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^ exp[-J2C{Vm,Vn)tmtn)e^Aj2tnVn \ T), (17) 

m,n n 

where C{Vm, Vn) are the two-body contact terms induced by the new 2-observables. Note 
that since these 2-observables are carried by the same 2-cycle 5, higher contact terms 
do not appear. The negative sign in front of the contact terms originates in the self- 
intersection number = —1 of the exceptional divisor. 

Some of these contact terms have been calculated explicitly. The simplest case of 
C{u2,U2) for G = SU{2) is due to Moore and Witten the result is written in terms 
of elliptic theta functions and Eisenstein series. This result is generalized to G = SU {N) 
by Marino and Moore 0. Their method is based on the so called RG (renormalization 
group) equation 

const. U2 (18) 



5 log A 

and the modular transformations 



(aiogA)2 (91ogA)2 91ogA9a/' ' ''d\og kda^ 

[iCr + Dr'],. ^_'[^^ (19) 



dlogkdaj dlogkdak 
under the symplectic transformations of cycles 

/ 



^3 ~^ CjkPk + Djkak, 



V 



A B 
C D 



ESp{2g,Z). (20) 



The second quantity above are the components of the directional vector V: 

const. . (21) 

o log Adaj 

Marino and Moore thus obtained the following formula: 

C{u2, U2) = const. jQ^^^- (22) 
By the results from the RG equation, this formula can be rewritten 

(Here the result is presented in a generalized form with Nf massless hypermultiplets). 
Losev et al. [§] derived a more general result: 

C^(«2,n.) = ^|P^(A:«.-a,|^). (24) 



(We have omitted an explicit form of various constants above, which are irrelevant in our 
subsequent analysis.) 



4 Multi-time tau function and contact terms 

Let us now proceed to our interpretation of the contact terms C{Vmi'Pn)- As we have 
noted above, the product of the last two factor in (p!oD persists to be a tau function of 
the Toda lattice hierarchy for any choice of the potential V of the 2-observable I{B). 
Meanwhile, the notion of the tau function can be readily generalized to the multi-time 
setting. Therefore, a natural conjecture is that right hand side of (|1^ will be a multi-time 
tau function of the Toda lattice hierarchy: 

T7,5(ti,t2,---) = exp(-5]C(P^,P„)t^t„)e^,5(E^n^" I (25) 

m,n n 

More precisely, this tau function should be an algebro-geometric tau function [|I^ 
determined by a suitable set of data (the Krichever data) on the Toda spectral curve (^) 
||TT| . (We here consider the case of Nf = only. Other cases can be treated in a similar 
way.) Such a tau function, just like the algebro-geometric tau functions of other integrable 
hierarchies takes the form 

r(ti, t2, ■■■) = e«(*^'*^'-)0(^ t„K + wo\T), (26) 

n 

where e*^'-*^'*^'" -' is a Gaussian factor (including a linear part), 

Q{ti,t2, ■■■) = ^Y1 Imntmtn + rjn + Tq, (27) 

m,n n 

<d{w I T) is the ordinary theta function without characteristic, and Wq is a constant 
vector. As we show later, qmn and Vn are determined by the the algebro-geometric data. 
r„, tq and Wq are arbitrary constants. Our conjectural tau function ( P^D can be rewritten 
into the above form using the relation 

e^^w I T) = e'"'<^'"'+'^>e(w + T^ + 5\T) (28) 

between the two theta functions. In particular, the arbitrary constants r„, Tq and Wq 
are determined by the half-characteristic and the period matrix, and qmn turn out to be 
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essentially the contact terms: 

(29) 

The above conjecture is supported by a modular transformation property of the tau 
function T^,<5(ti, ^2, ■ ' ') under symplectic transformations of cycles (pO]). Namely, as we 
shall show below, the tau function turns out to possess the same t-INDEPENDENT mod- 
ular property as the original t-dependent factors in (|10|) ; this implies that the corrected 
0-observable factor (^) persists to be modular invariant after modifying the last two 
factors as in (|1^. This is strong evidence (though not a proof) of the correctness of the 
conjecture. 

Let us examine modular transformations of the tau function T^,5(ti, ^2? ■ ■ ■)• Fortu- 
nately, this kind of problems have been already studied in the context of free fermion 
systems on Riemann surfaces |]T3|. (This also suggests a possible link of four- dimensional 
super symmetric gauge theories with two-dimensional conformal field theories.) Remark- 
ably, the modular transformation property of the tau function is quite universal, i.e., does 
not depend on the detail of the Riemann surface or the algebro-geometric data (even nor 
the integrable system itself!). In the following, therefore, we consider the tau function 
^7,5(^1,^2, ■ ■ ■) for an arbitrary compact Riemann surface S of genus g. 

The fundamental algebro-geometric data in the case of the Toda lattice hierarchy [|1^ 
are comprised of two marked points P^, local complex coordinates z± at each of those 
points normalized as z±{P^) = 0, and the symplectic basis of cycles aj, (3j {j = 1, - ■ ■ ,g). 
Furthermore, a set of polynomials f^{z^^) {n = 1,2,---) have to be given in order to 
select the "directions" of the t„'s in the total space of the standard time variables 
(n = 1, 2, ■ ■ ■) of the Toda lattice hierarchy 0. Each pair f:^{z^^) and f~{z^^) determines 
the direction of the n-th time (or, equivalently, a Hamiltonian Vn in the sense already 
mentioned) . 

Given these data, the following meromorphic differentials dfln = 1, 2, ■ ■ ■) are 
uniquely determined: 

• dQn is holomorphic everywhere except at P^. In a neighborhood of P^, respectively, 

d^n = df^{zg^) + holomorphic. (30) 
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• The a-periods all vanish, 

/ rffi„ = {j = l,---,g). (31) 

These meromorphic differentials determine qmn and Vn as follows. The components of 
the vector Vn are given by 



Vjn = — f d^n. (32) 
2711 Jf3j 

By Riemann's bilinear relation, this can be rewritten 

where the integrals on the right hand side are a contour integral turning once around 
respectively, q^n is given by replacing dojj by dVLn- 

Qmn = f^{z^^)dQrn ~ ^ fni^-^)d^m- (34) 

These somewhat complicated expressions are a linear combination of more familiar ex- 
pressions of the g's and V^s for the standard time variables [10|. 



The modular transformation properties of these quantities can be derived by straight- 
forward calculations. Under the symplectic transformation of cycles (^), the holomorphic 
and meromorphic differentials transform as 

dujj [{CT + D)~^]kjduk, 

dQn ^ dnn-[{CT + D)-'C]ke(^f^dn^^duJk. (35) 
Accordingly, Vn and Qmn transform as 

Vn ^ \CT + Dr'Vn, 
Qmn ^ qn.n - 2m'Vn,{Cr + D)-^CVn. (36) 

The final piece of the ring is the following modular transformation formula of theta 



functions [14 



e^,s{\CT + D)-^w I {AT + B){CT + D)-^] 
= edet{CT + DY^^expi-K < w,\CT + D)-^w >)Qy^s'iw \ T). (37) 
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Here e is an 8**^ root of unity, = 1, and (7', 5') is a new half-characteristic, both 
determined by the symplectic matrix. Combining this formula with the above modular 
transformations of qmn and Vn lead to the conclusion that the tau function r^,5(ti, ^2, ■ ' ') 
has a t-independent modular transformation property: 

T^,sih,t2,---)^edetiCT + DY'\,^s'{tiM,- ■ ■)■ (38) 

Understanding that the half-characteristic is also transformed because of its physical 
origin, one can thus confirm that the tau function in our conjecture possess a desirable 
modular transformation property. 



5 Perspectives from Whitham equations and prepo- 
tent ial 

Let us reconsider the meaning of Qmni Vjn and Vn in the language of Whitham deformations 
of the integrable hierarchy. The most suggestive in this respect is the formula (p2|) of 
C{u2,U2) as a second derivative of the prepotential. The relation between RG equations 
and Whitham equations tells us that, very roughly speaking, the energy scale A can 
be identified with the first Whitham time variable Ti. In view of this fact, a natural 
conjecture is that the contact terms of higher observables can be written 

C(P™,P„) = const. (39) 



This conjecture indeed turns out to fit into the general framework of Whitham equations 
proposed in our previous work as we show below. 

The Whitham equations in the present situation take the form 

^T„C^>S'|2=const. = dVLn^ dajdS\z=const. = dUj, (40) 

where dS is given by 

N-l 

dS = '^TndQn + o-jd^j^ (41) 

n j=l 

and "|2=const.", also here, means differentiating while keeping z constant. These equations 
give deformations of the Coulomb moduli Uj = Uj{a,T) as the Whitham time variables 
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T = (Ti,T2,- ■ ■) vary from a point of departure (e.g., a point T = Tsw where dS coincides 
with the meromorphic differential dSsw)- One can redefine the prepotential JF = jF(a, T), 
now as a function of a and T, by the equations 

f^{z-')dS- i J-(zZ')dS (42) 



dT Jp+ + /p 
and 

dS. (43) 

Recall that f^{z^^) are the polynomials giving the singular part of dVLn at P^. The com- 
patibility (integrability) of the above equations for JF is again a consequence of Riemann's 
bilinear relation. 

Now differentiate the right hand side of the defining equation of dT jdTn against aj 
and T^. By the Whitham equations, the outcome is nothing but the right hand side of 
(p3|) and (p^). Thus Vj„ and qmn-, which are now functions of a and T, can be expressed 
as second derivatives of the prepotential: 

1 d^T 1 d^T 

" 2^%d^n' " 2mdf^- ^^^^ 

Similarly, the matrix elements of the period matrix can be written 

'^^ da j dak' ^^^^ 
Since Vn should be the gradient vector of the potential Vn in the a-space, we conclude 
that the potential Vn can be written 

Thus all relevant quantities turn out to be written as derivatives of the redefined prepoten- 
tial. It is remarkable that the last equation resembles the hypothetical Hamilton- Jacobi 
equation of Losev et al. In our case, however, Vn is a function of both the Coulomb 
moduli AND the Whitham time variables, the latter enter from dS through the integral 
formula, and we do not know how to convert it into a function of aj and af = dT/daj 
like the Hamiltonians of Losev et al. 

The redefined prepotential is a purely theoretical backbone, and not very suited for 
explicit calculations. Integral formulae, such as (0), ( ^4|) and an integral representation 
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of Vn derived from (|4^) , are more convenient. These integral formulae remain valid even 
if the Whitham time variables are returned to the "departure time" T = Tsw where dS 
is equal to dSsw] everything can be thereby calculated in terms of the original geometric 
data on the spectral curve. Calculations of contact terms are thus eventually reduced to 
residue calculus. This is enough for understanding the four-dimensional problem. 

Nevertheless, special families of Whitham deformations can possess some significant 
implications. Gorsky et al. indeed presented such an example in the case of G = 
SU (N) and Nf = 0. As they remarked, their Whitham deformations exhibit a remarkable 
similarity with two-dimensional topological Landau- Ginzburg models [16]. The construc- 
tion starts from the meromorphic differentials 

dQn = dlogz, (47) 



where (■ ■ ■)+ means the polynomial part of the Laurent expansion of P(x)"''^ at x = oo. 
dCli is nothing but the differential dS^w- Now introduce a set of time variables T„ and 
consider the differential 

oo 

dS=Y^ TA. (48) 

n=l 

Because of a reason (see below), we have to distinguish between these time variables T„ 
and the previous ones T„. The period integrals 

aj = <f dS. (49) 

define a function of the Coulomb moduli u = {u2, ■ ■ ■ ,un) and the time variables Ti, T2, ■ ■ ■. 
One can prove, by a standard method in Seiberg-Witten geometry, that the Jacobian 
matrix det{daj/duk) does not vanish in a neighborhood of, say, the "Seiberg-Witten 
point" Tn = 6n,i- Therefore the above relation can be solved for the Coulomb moduli 
as Uj = Uj{a,Ti,T2, ■ ■ ■). This gives a deformation family of the spectral curve S. Now 
modify the meromorphic differential dCln into 

N-l 

dVln = dCln ~ X! ( r dVlnjdujj, (50) 

and write dS in the form 

00 ^ N-l 

dS = TndCln + ^ QjdUj. (51) 

n=l j=l 
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One can then derive, by the method of Itoyama and Morozov 
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, the Whitham equations 



df^dS\z=const. — dVtn, dajdS\z=const. — dUj. (52) 

This Whitham deformation family is shghtly distinct from those that we have consid- 
ered. A natural choice of local coordinates Zj- at is the following: 

z+ = z-i/^, z_ = z^"". (53) 

The polynomials /^(^^^) giving the singular part of d^n at can be written 

n 

UzZ') = --A"zZ" + .... (54) 
n 

The difference lies in the tail part "■ ■ ■" . This part vanishes for n < 2N, but remains for 
n > 2N, and the coefficients of this part are NOT a numerical constant but a polynomial 
of the moduli Wj's. This is the reason that we changed the notation of the time variables. 
This somewhat strange situation forces us a careful treatment of the prepotential. For 
instance, when differentiating contour integrals like those in ( ^2]) against and a^., we 
use the fact that the part of f^{z^^) may be considered constant; this is not permitted 
in the above example. 

Extending this example to other gauge groups, such as S0{2N), is an interesting 
problem. We shall consider this issue elsewhere. 



6 Conclusion 

Inspired by the work of Marino and Moore [0, we have proposed an extension of the 
blowup formula with more than one 2-observables of the form /„(-B) supported on the 
exceptional divisor B. Our strategy is simply to replace the single-time tau function of the 
Toda lattice hierarchy (in the sense of Marino and Moore) by a multi-time tau function. 
The time variables t„ are interpreted as the coupling constants for the insertion of In{B). 
The tau function is an algebro-geometric tau function comprised of a Gaussian factor 
gQ(ii,i2,.--) QT^d a theta function 07,5(I]n '^nKi I ^)- The coefficients of the Gaussian part 
are identified to be the contact terms of the 2-observables. We have partly confirmed the 
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validity of our proposal by showing that the tau function possesses a desirable modular 
transformation property. 

This proposal has been further examined in the language of the Whitham equations 
that underlie the integrable hierarchy. We have shown that the contact terms, as well as 
other relevant quantities, are written in the form of derivatives of the prepotential with 
respect to the Whitham time variables T„. This also clearly explains why the RG equation 
takes place in the description of contact terms of the quadratic Casimir U2- 

These observations should be further checked in a field-theoretic language. The su- 
perfield formalism of Losev et al. 0| will provide a suitable framework for this purpose. 

It is remarkable that the integrable hierarchy and the Whitham hierarchy are both 
linked with two-dimensional field theories. The integrable hierarchy (the Toda lattice 
hierarchy in the case of G = SU{N) and Nf = 0) is related to massless free fermions 
on the spectral curve [jl8|. (This is also an implicit message from the work of Gorsky et 
al. They used the Szego kernels, which are correlation functions of free fermions.) 
The modular transformation property of the tau function is physically a consequence of 
conformal invariance of the two-dimensional massless free fermion theory. Meanwhile, 
as recent attempts at an analogue of the WDVV equations [1^ suggest, the Whitham 
hierarchy are related to two-dimensional topological CFT's, in particular, topological 
Landau- Ginzburg models. 

These two-dimensional structures deserve to be studied in more detail. Of particular 
interest will be to examine, from our point of view, the mirror-like structure pointed out 



by Ito and Yang [IS]. 



Acknowledgements 

The author is grateful to Toshio Nakatsu for discussions. He also thanks Masahiko Saito 
for invitation to his workshop held in Kinosaki, a beautiful hot spring town, where this 
work was started. This work is partly supported by the Grant-in-Aid for Scientific Re- 
searches (No. 09304002), the Ministry of Education, Science and Culture, Japan. 



References 



14 



[1] G. Moore and E. Witten, Integration over the w-plane in Donaldson Theory, Piep- 
th/9700193. 



[2] M. Marifio and G. Moore, Integrating over the Coulomb branch in N=2 gauge theory, 
Nucl. Phys. Proc. Suppl. 68 (1998), 336-347; The Donaldson- Witten function for 
gauge groups of rank greater than one, Commun. Math. Phys. 199 (1998), 25-69. 

[3] A. Losev, N. Nekrasov and S. Shatashvili, Issues in topological gauge theory, Nucl. 
Phys. B534 (1998), 549-611; Testing Seiberg- Witten solution, |hep-th/98U 10611 . 



[4] E. Witten, Topological quantum field theory, Commun. Math. Phys. 117 (1988), 
353-386; Supersymmetric Yang-Mills theory on a four-manifold, J. Math. Phys. 35 
(1994), 5101-5135. 

[5] N. Seiberg and E. Witten, Electric-magnetic duality, monopole condensation, and 
confinement in N = 2 supersymmetric Yang-Mills theory, Nucl. Phys. B426 (1994), 
19-52; Monopoles, duality and chiral symmetry breaking in N = 2 supersymmetric 
QCD, Nucl. Phys. B431 (1994), 484-550. 

[6] A. Gorsky, I. Krichever, A. Marshakov, A. Mironov and A. Morozov, Integrability 
and Seiberg- Witten exact solution, Phys. Lett. 335B (1995), 466-474. 
E. Martinec and N.P. Warner, Integrable systems and supersymmetric gauge theory, 
Nucl. Phys. B459 (1996), 97-112. 

T. Nakatsu and K. Takasaki, Whitham-Toda hierarchy and N = 2 supersymmetric 
Yang-Mills theory. Mod. Phys. Lett. All (1996), 157-161. 

R. Donagi and E. Witten, Supersymmetric Yang-Mills theory and integrable systems, 
Nucl. Phys. B460 (1996), 299-344. 

A. Gorsky and A. Marshakov, Towards effective topological gauge theories on spec- 
tral curves, Phys. Lett. B375 (1996), 127. 

H. Itoyama and A. Morozov, Integrability and Seiberg- Witten theory, Nucl. Phys. 
B477 (1996), 855-877; Prepotential and the Seiberg- Witten theory, Nucl. Phys. 
B491 (1997), 529-573. 

C. Ahn and S. Nam, Integrable structure in supersymmetric gauge theories with 
massive hypermultiplets, Phys. Lett. B387 (1996), 304-309. 



15 



I.M. Krichever and D.H. Phong, On the integrable geometry of sohton equations and 
N = 2 supersymmetric gauge theories, J. Diff. Geom. 45 (1997), 349-389. 

[7] A. Klemm, W. Lerche, S. Theisen and S. Yankielowicz, Simple singularities and N = 
2 supersymmetric Yang-Mills theory, Phys. Lett. 344B (1995), 169-175. 
P. Argyres and A. Faraggi, The vacuum structure and spectrum of N = 2 supersym- 
metric SU(n) gauge theory, Phys. Rev. Lett. 73 (1995), 3931-3934. 

[8] K. Ueno and K. Takasaki, Toda lattice hierarchy. Adv. Stud. Pure Math. 4 (1984), 
1-95. 

T. Takebe, Toda lattice hierarchy and conservation laws, Commun. Math. Phys. 129 
(1990), 281-318. 

[9] G. Bonelli and M. Matone, Nonperturbative renormalization group equation and 
beta function in N = 2 SUSY Yang-Mills, Phys. Rev. Lett. 76 (1996), 4107-4110. 
E. D'Hoker, LM. Krichever and D.H. Phong, The renormalization group equations 
in N = 2 supersymmetric gauge theories, Nucl. Phys. B494 (1997), 89-104. 
E. D'hoker and D.H. Phong, Calogero-Moser systems in SU(N) Seiberg-Witten the- 
ory, Nucl. Phys. B513 (1998), 405-444. 

A. Gorsky, A. Marshakov, A. Mironov and A. Morozov, RG equations from Whitham 
hierarchy, Nucl. Phys. 527 (1998), 690-716. 

[10] T. Nakatsu and K. Takasaki, cited above 0. 

[11] LM. Krichever, Algebraic curves and non-linear difference equations, Russ. Math. 
Surveys 33:4 (1978), 215-216; Theta functions and non-linear equations, Russ. Math. 
Surveys 36:2 (1981), 11-92 (Appendix). 

[12] M. Jimbo and T. Miwa, Monodromy preserving deformations of linear ordinary dif- 
ferential equations with rational coefficients II, Physica 2D (1981), 407-448; ditto 111, 
Physica 4D (1981), 26-46. 

E. Date, M. Jimbo, M. Kashiwara and T. Miwa, Transformation theory for soliton 
equations HI, J. Phys. Soc. Japan 50 (1982), 3806-3812; ditto IV, J. Phys. Soc. Japan 
50 (1982), 3813-3818; ditto V, Physica 4D (1982), 343-365; ditto VI, Publ. RIMS., 
Kyoto Univ., 18 (1982), 1077-1110. 

16 



T. Shiota, Characterization of jacobian varieties in terms of soliton equations, Invent. 
Math. 83 (1986), 333-382. 

[13] L. Alvarez- Gaume, G. Moore and C. Vafa, Theta functions, modular invariance and 
strings, Commun. Math. Phys. 106 (1986), 1-40. 

N. Kawamoto, Y. Namikawa, A. Tsuchiya and Y. Yamada, Geometric realization 
of conformal field theory on Riemann surfaces, Commun. Math. Phys. 116 (1988), 
247-308. 

[14] D. Mumford, Tata lectures on theta (Birkhauser, Boston, 1983). 

[15] A. Gorsky, A. Marshakov, A. Mironov and A. Morozov, cited above [Q. 

[16] Dijkgraaf, R., Verlinde, E., and Verlinde, H., Topological strings in c? < 1, Nucl. 
Phys. B352 (1991), 59-86. 

[17] H. Itoyama and A. Morozov, cited above 

[18] L. Alvarez-Gaume, C. Gomez and C. Raina, Loop groups, grassmannians and string 
theory, Phys. Lett. B190 (1987), 55. 

N. Ishibashi, Y. Matsuo and H. Ooguri, Soliton equations and free fermions on Rie- 
mann surfaces. Mod. Phys. Lett. A2 (1987), 119. 

C. Vafa, Operator formulation on Riemann surfaces, Phys. Lett. B190 (1987), 47. 

[19] G. Bonelli and M. Matone, Nonperturbative relations in N = 2 SUSY Yang-Mills 
and WDVV equation, Phys. Rev. Lett. 77 (1996), 4712-4715. 

A. Marshakov, A. Mironov and A. Morozov, WDVV-like equations in N = 2 SUSY 
Yang-Mills Theory, Phys. Lett. B389 (1996), 43-52; WDVV equations from algebra 
of forms. Mod. Phys. Lett. A12 (1997), 773-788. 

K. Ito and S.-K. Yang, The WDVV equations in N = 2 supersymmetric Yang-Mills 
theory, |hep-th/9803T26 . 



17 



